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Abstract 

We compute the time evolution of the mutual information in out of equilibrium quantum systems 
whose gravity duals are Vaidya spacetimes in three and four dimensions, which describe the 
formation of a black hole through the collapse of null dust. We find the holographic mutual 
information to be non monotonic in time and always monogamous in the ranges explored. We 
also find that there is a region in the configuration space where it vanishes at all times. We show 
that the null energy condition is a necessary condition for both the strong subadditivity of the 
holographic entanglement entropy and the monogamy of the holographic mutual information. 
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Introduction 



Entanglement entropy is a measure of the quantum entanglement in systems with many de- 
grees of freedom which has been object of intense investigation in condensed matter, quantum 
information and quantum gravity. 

Given a system whose total Hilbert space H can be written as a direct product H = Ha®Hb 
and whose state is characterized by the density matrix p, one can define the reduced density 
matrix pa = Tr#/? by tracing over the degrees of freedom of B. Then, the entanglement entropy 
Sa is the corresponding Von Neumann entropy Sa = —T^a(pa log pa)- The the situation mostly 
studied in the literature is when A is given by a spatial region and B is its complement. In 
this case the entanglement entropy is also called geometric entropy and Sa behaves according 
to an area law: in d > 1 spatial dimensions we have Sa oc Area(cM)/e d_1 + . . . , where e is the 
UV cutoff of the theory and the dots represent terms of higher order in e pQ . In two spacetime 
dimensions the divergence is logarithmic and more quantitative analysis has been performed for 
conformal field theories, where the symmetry provides powerful computational techniques. For 
these theories, when A is given by an interval of length I in an infinite line and the temperature 
is zero, we have the famous result Sa = (c/3) log(£/e) + c' l5 where c is the central charge and c' x 
is a constant [21 El E] (see [5] for a recent review). The most useful method to get Sa is the so 
called replica trick, which consists in computing the Renyi entropies = (1 — n) _1 logTr/?^ 
for integer n. The entanglement entropy is then obtained as Sa = — d n TrpW n= \. 

For quantum field theories with a holographic dual, the prescription to compute Sa in the 
boundary theory through a bulk computation has been conjectured in [6l [7] for static back- 
grounds and then generalized to stationary and time dependent geometries in [8] (see [9] for a 
review). The proposal of [6 , 7J satisfies many properties and also important inequalities of the en- 
tanglement entropy (the simplest of them is the strong subadditivity) |10l lilt [T2] . Nevertheless, 
a proof for this formula is not known |13| W$\ . 

Interesting insights in the structure of entanglement can be obtained by considering two 
disjoint regions, namely A = A\ U A2 with A\ n A2 = [US EH El [IB]. In this case the 
most interesting quantity to study is the mutual information I(A\,A2) = 5a 1 + Sa 2 ~ 5a 1 uA 2 > 
which has the nice property that the leading divergence due to the area law cancels. The Renyi 
mutual information has been analytically computed for some simple two dimensional conformal 
field theories like the free compactified boson |17] and the Ising model [18J, and it has been 
shown that it encodes all the data of the theory (all the conformal dimensions of the primaries 
and their correlation functions) in contrast with the entanglement entropy of a single interval 
which contains only the central charge. Unfortunately, the analytic continuation that allows to 
obtain I(A\, A2) is not known in general but only in some limiting regimes. Detailed studies of 
I(Ai, A2) for spin chain models have also been done [HI EUJ EH [22] . 

The holographic formula of [HE] for static backgrounds has been applied for disjoint regions 
|23l HHl [23] and a qualitative disagreement with respect to the expectations from the simple two 
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dimensional conformal field theories mentioned above has been found. Indeed the holographic 
mutual information displays a continuos transition with discontinuous first derivative from zero 
to positive values as the two regions get closer. This could be explained through the fact that 
the holographic formula holds for large c and it should be corrected in order to recover the 
results obtained for the models whose central charges are of the order of the unity. Nevertheless 
the holographic formula of E] and its generalization for time dependent backgrounds [S] are 
believed to be correct for large c and we will employ this prescription in our analysis. 

The entanglement entropy is a very important quantity to study in order to understand the 
physics out of equilibrium and the processes of thermalization. In particular, one is interested 
in the unitary time evolution of the entanglement entropy when the system starts from a state 
which is not an eigenstate of the Hamiltonian of the system. This occurs for instance when a 
system is prepared in an eigenstate of the corresponding Hamiltonian and suddenly a tunable 
parameter of the Hamiltonian (e.g. the magnetic field) is changed (global quench). Then the 
system evolves unitarily accordingly with the new Hamiltonian starting from a state which is 
not one of its eigenstates. Another interesting situation is when the system is prepared in 
the ground state of two decoupled parts which are joined together at t = and then the whole 
system unitarily evolves through the translationally invariant Hamiltonian (local quench). These 
situations have been studied for two dimensional systems by employing conformal field theory 
techniques and spin chains models [251 ESI EH ES EH] • 

Thermalization processes have been also widely considered from the holographic point of 
view (see [30J for a review). For instance, the holographic counterpart of the unitary evolution 
of a system towards a stationary state is a gravitational collapse whose final state is a black hole. 
In this process both the initial and the final states are thermal and the holographic entanglement 
entropy can be computed through the prescription of [BJ. 

The Vaidya metrics (see e.g. |31] ) are simple backgrounds realizing this setup. They are so- 
lutions of the Einstein equations with negative cosmological constant and a non trivial energy 
momentum tensor containing a mass function constrained by the null energy condition which 
describe the formation of a black hole through the collapse of a shell of null dust. The null 
geodesies in these geometries have been studied in j32j but here we are interested into the space- 
like ones, which occur in the computation of the holographic entanglement entropy [U [33]. I34j . 
The Vaidya metrics are simplifications of more general models considered e.g. in [35] (tensionful 
shell) and [36] . Other holographic thermalization setup have been also studied [37] . The study of 
the holographic entanglement entropy in Vaidya backgrounds is usually numerical but recently 
an analytic computation has been done for the limiting regime of thin shell in three spacetime 
dimensions |38[ 139] . We will largely employ this result in our analysis because it allows to ex- 
plore a larger range of parameters. The holographic analysis of the two point functions in these 
dynamical geometries has been carried out in |38[ [39l 130] . 
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In this paper we study the holographic entanglement entropy in the three and four di- 
mensional Vaidya backgrounds. In section [I] we introduce the metrics and the corresponding 
holographic entanglement entropy for a single region in the boundary theory. In section [2] we 
study the holographic mutual information and its transition curves in the configuration space. 
In section [3] we explore the relation between the null energy condition for the Vaidya metrics 
and the strong subadditivity for the holographic entanglement entropy. In section [4] we extend 
the analysis performed in the previous sections to the holographic tripartite information in order 
to verify the monogamy of the holographic mutual information in Vaidya spacetimes and study 
this property is influenced by the null energy condition. 

Note added. While we were completing the writing of this paper, [3T] appeared and it has a 
substantial overlap with our results. 



1 Holographic entanglement entropy for Vaidya geometries 



In this section we introduce the Vaidya metrics in d + 1 dimensions (subsection and we 
describe some known results about the holographic entanglement entropy in these backgrounds 
(subsection 1.2). In subsection 1.3 we focus on the three dimensional case because it is the 
simplest to study and some analytical results have recently been found |38^ [3"5] which will be 
widely employed in the remaining part of the paper. 



1.1 Vaidya metrics 

The d + 1 dimensional Vaidya metrics in the Poincare coordinates read 



ds 2 



2 



Z 



m(v)z d )dv 2 - 2dzdv + dx 1 



(1.1) 



(we set &itG^ +1 ^ = 1) where x = {x\, . . . , Xd-i} are the spatial boundary coordinates. The 
Ricci scalar of (1.1) is R = — (d + l)d/l 2 . The metric (1.1) is a solution of the Einstein equations 
in presence of matter 



Gnu + Ao 



T, 



A 



d(d + l) 
2P~ 



where the energy momentum tensor has only one non vanishing component 



d 



- z^dvmiv) 



r X2) 



[1.3) 



The metric ( 1.1 ) describes the formation of a black hole through the collapse of null dust, which 
is characterized by the just introduced. It is important to observe that the translational 
invariance along the directions contained in x is preserved at any time. This is a key feature of 
the setup characterizing global quenches in the boundary theory. 

In processes of gravitational collapse, it is not yet understood how to characterize the formation 
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of a black hole through local time evolution. To this purpose, some generalizations of the 
concept of horizon have been proposed, and, in our backgrounds we can distinguish between two 
horizons, the event and the apparent horizon [42 } 143 1 EH]. The apparent horizon is the boundary 



of the trapped surfaces associated to a given foliation. For the metrics (1.1) a foliation which 



preserves the translation invariance in the directions of x is given by v = const and z = const. 



The location of the apparent horizon of (1.1) reads [51 133] 

1 



m(v 



Instead, the event horizon is given by 



dz e 
dv 



m{v)z\ 
2 



(1.4) 



(1.5) 



When the mass profile m{v) is constant m(v) = M, the metric (1.1) describes the geometry of 



the Schwarzschild black hole with planar horizon. This can be clearly seen through the following 
change of coordinates 



t + p(z) 



p'(z) 



1 



1 - Mz° 



which allows to write (1.1) as 



ds 2 



z 



dz' 



[\ - Mz d )dt 2 + — , + dx 

K ' 1 - Mz d 



(1.6) 



(1.7) 



i.e. the usual form for a Schwarzschild black hole of mass M in the Poincare coordinates. For 
I = 1 the Hawking temperature of this black hole is given by T H = dM l / d /{ATt). It is also 



straightforward to check that for m(v) = identically the metric (1.1) describes AdSd+i in the 
Poincare coordinates 



ds 1 



I 2 



dt 2 + dz 2 + dx 2 ) 



t 



V + z . 



(1.8) 



This tells us that the Vaidya metric (1.1) is asymptotically AdSd+i- 
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1. The dashed 

curve is a step function and corresponds to the the limiting regime of thin shell a v — > 0. 



Figure 1: The function (1.9) for different values of the thickness a v and M 
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The formation of the Schwarzschild black hole (1.7) from AdSd+i is described by considering a 



function m(v) which interpolates between and the finite value M > in an strictly increasing 
way. The profile for m(v) is usually chosen to be 

. , „ , 1 + t&nh.(v/a v ) 
m(v) = M v ' 1 . (1.9) 



Given this m(v), the metric ( 1.1 ) describes a background which evolves from pure planar AdSd+i 
at early times to Schwarzchild black brane with mass M at late times because of the infalling 
shell of null dust. The parameter a v determines the transition between these two regimes since it 



parameterizes the thickness of the shell which falls along v = 0. The mass function ( 1.9 ) is shown 
in the figure [TJ In the limiting case of a v — > 0, the mass profile m(v) becomes a step function 
M6(v) and the infalling shell describes a shock wave. This limit is very interesting because it 
captures the essential physics of the problem and one can hope to find analytic solutions for the 
quantities considered. For the holographic entanglement entropy in three spacetime dimensions 
this has been done in [38, 39J and we will largely employ this result. 



An important inequality to impose on the energy momentum tensor in (1.2) to guarantee 
the positivity of the energy density is the null energy condition, namely T^N^N" ^ for any 
null vector |42|. H3] . This inequality has been employed to study the c theorems from the 
holographic point of view [45(. |4T>] and their generalizations in presence of boundaries |47l 0H1 



50] . For the energy momentum tensor (1.3), imposing the null energy condition means 



d v m(v) > (1.10) 



which is clearly satisfied by the profile (1.9). In the sections [3] and [4] we consider mass profiles 



violating this condition and the effect of this violation on the holographic entanglement entropy. 
1.2 Holographic entanglement entropy 



For static spacetimes like AdS^+i (1.8) and the Schwarzschild black hole (1.7), the prescription 
to obtain the entanglement entropy Sa = — Tr(pA log pa) in the boundary theory through a 
holographic computation in the bulk has been proposed in (6j [7]. It reads 

4G^ +1) 

where G^ +1 ^ is the Newton constant for the d + 1 dimensional bulk and ja is defined as the 
minimal surface among the spatial ones which extend into the bulk and share the boundary 
with A, i.e. d^A = 9A (■ja is homologous to A). Thus, ja is a codimension two surface living in 
the t = const slice and has minimal area. Since 7^4 lives in an asymptotically AdSd+i space and 
it reaches its boundary, placed at z = in the Poincare coordinates, Area(7^) is infinite and 
therefore it must be regularized by introducing a small cutoff e > in the holographic direction, 
namely the restriction z > e is imposed during the bulk computation. The series expansion of 
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5,4 in e depends on d and one the first checks of was that this expansion reproduces in the 

holographic way the leading UV divergence of the entanglement entropy which was computed 
in field theory through other methods. In particular, for d > 2, the leading divergence of Sa 
is proportional through a non universal coefficient to Area(<9 A) / 'e d ~ 2 [I] (this is the so called 
area law). For d = 2, when the boundary CFT is two dimensional and the spatial region 
A is a segment of length £ (thus dA is made by two points), the entanglement entropy Sa 
diverges logarithmically with a universal coefficient given by the central charge of the theory 
[21 El Sj. Besides these fundamental checks, it has been shown that the holographic proposal 
satisfies the strong subadditivity condition and other more complicated inequalities 
characterizing the entanglement entropy [12] , which will be discussed in sections [3] and |4j 

The proposal ( |1.11[ ) for static backgrounds has been generalized to time dependent geome- 
tries in [8]. In these cases Sa is still given by ( |1.11[ ) but with ja obtained as an extremal surface, 
namely as the saddle point of the proper area functional. This proposal is covariantly well de- 
fined and reduces to the previous one when the spacetime is static. Let us explain it in details 



for the Vaidya d + 1 dimensional spacetimes (1.1) 



We consider as d — 1 dimensional spatial region A in the boundary theory a rectangular strip 
parameterized by x\ S (—£/2, £/2) and X2, ■ ■ ■ Xd-i £ (0, £±) at fixed value of the boundary time 
coordinate t. This choice is less symmetric than the case of A given by a circular shape but a 
crucial simplification occurs in this case. 

According to the proposal of [8], the holographic entanglement for this spatial region A is given 
by the area of the extremal surface 7A whose profile is most conveniently specified by v = v(x\) 
and z = z(x\) (we are assuming that 7^ is translationally invariant in the other boundary 
coordinates X2, . . . Xd-i parameterizing A) with the following boundary conditions 

v{-£/2) = v{£/2) =t z(-£/2) = z{£/2) = . (1.12) 

These boundary conditions impose that the boundary of ja coincides with the boundary of A 
along the boundary temporal evolution. Since x\ is the relevant independent variable, we will 
denote it by x in the following. The area of the such spacelike surface (thus the determinant of 
the induced metric under the square root must be taken with the positive sign) is given by the 
following functional 

Area( 7A ) = A d = l^ 1 2£ d ~ 2 [" J I - [l - m(v)z d ] {v') 2 - 2v'z' dx (1.13) 

Jo 2 v 

where ' = d/dx and we have chosen the origin along the x direction in order to employ the fact 
that the functions v(x) and z(x) are even. This consideration about parity determines the factor 



2 and the integration extrema in (1.13). 



The area functional ( 1.13 ) is the one we have to extremize in order to get the the codimension two 
surface 7^ allowing to compute the holographic entanglement entropy for the time dependent 
Vaidya spaces through ( |1.11[ ) [SJ. In other words, ja is a solution of the two equations of motion 
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of (1.13). Since the integrand in (1.13) does not contain x explicitly, we have the following 

[1 - m{v)z d ](v'f 



conservation equation 

■j'z' (1.14) 

where z* = z(0) is the maximum value of z(x), which is characterized by z' = v' = 0. The 
occurrence of this conservation law is the simplification characterizing for the rectangular shape 
of A with respect to the circular one mentioned above. The two equations of motion obtained 



by minimizing the functional (1.13) read 

[1 - m(v)z d ]v" + z" 
„ d-2 



di ii i // d v m(v) d r 2 j , s d-l i i 
" — z [v ) — dm[y)z zv 



Z V 



m{v)z d {v'f + {d - 1) [ {v'f + 2v'z' - 1 ] 



. 



(1.15) 
(1.16) 



By taking the derivative w.r.t. x of the conservation equation (1.14) and using one of these two 



equations of motion, one obtains the other one. Thus, it is sufficient to consider only (1.14) and 



e.g. (1.16) to find v(x) and z(x). 
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Figure 2: Finite term of the holographic entanglement entropy for the Vaidya metric in three 



dimensions (see (1.11) and (1.18) up to a factor) with m(v) given by (1.9) [U|33]. On the left, 
-^2,reg(^)£) for different values^ of boundary time t, increasing from the red curve to the blue 
one, at a fixed value of a v . On the right, L2,reg(^> o-v) at a fixed value of t and different values of 
the thickness a v , which increases going from the blue curve to the red one. In both the plots, 



the black curves correspond to the limiting regimes of AdS% (bottom curve, from (1.24)) and of 



the BTZ blak hole (top curve, from (1.24)), while the dashed curves represent the corresponding 
curves for the thin shell limit a v = 0. 

Given a solution {v (x), z{x)) of the equations of motion, its area can be computed by evaluating 
the integral ( 1.13[ ) on it. Using (1.14), this area can be written as follows 



A d = I 



d—l 2£<i— 2 



z 2(d-l) 



dx . 



(1.17) 



'In all the plots of this paper the numerical values of the parameter characterizing the different curves are 
evenly spaced. 
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Figure 3: Finite term of the holographic entanglement entropy for the Vaidya metric in four 
dimensions (proportional to L3 ireg (£, t)) with m(v) given by ( 1.9 ), fixed a v and different boundary 



times t, increasing from the red curve to the blue one. The black curves correspond to the 



limiting regimes of AdS4 (bottom curve from (1.20)) and of the Schwarzschild black hole in four 
dimensions (top curve). 



As discussed above, this integral is divergent because the spacetime we are dealing with, being 
asymptotically AdSd+i, is non compact and the spatial surface we are considering reaches 



its boundary (see the boundary conditions (1.12)). 

The divergence of Ad can be obtained by studying the same problem in AdSd+i in the standard 
way [7j- Subtracting this divergence we obtain the finite term of the area which is the main 
quantity we are interested in. For d > 2 we have 



A d . 



21 



d — 1 £d — 2 



lim 

r;^0+ 



, -77 „d— 1 



z 2(d-l) 



dx 



1 



(d-2)e 



d-2 



l d - l Z d ,- 2 L d . 



reg 



;i.i8) 



where e = z(£/2 — rf) is the UV cutoff in the boundary theory. Notice that in this case we need 
to subtract just one diverging term to regularize A d - This is a feature of the strip; indeed when 
the region A is a circle there are more terms to subtract to make the area finite [7]. 



In order to find the solution of (1.14) and (1.16) satisfying the boundary conditions (1.12), 



first we exploit the reflection symmetry about x = and solve the Cauchy problem whose initial 
conditions are given by 

z(0) = v(0) = v* . (1.19) 

Then, we shoot in the variables z*, v* to impose (1.12F It turns out that the points at which 
the solution reaches the boundary become increasingly sensitive to the initial conditions and to 



"This approach is different from the one used in [33| . where shooting takes place from the boundary. We 
favored shooting from x — because it is a regular point of the solution (v(x), z(x)). 
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the accuracy of the integration as either d or £ increase. As a consequence, it becomes more 



and more difficult to impose (1.12). This technical difficulty limits the range of parameters we 
are able to explore for d > 2. Once the solution is found, the implementation of the numerical 
integration and limit in ( |1.18[ ) is quite straightforward. 

In the figures [2] and [3] we show respectively ^2,reg (three dimensional Vaidya) and ^3,reg 
(four dimensional Vaidya) as functions of £ and for different values of the two other important 
parameters involved in our problem: the boundary time t and the thickness of the shell a v . 
The three dimensional case deserves a separated discussion (see subsection 1.3) because it is 
the simplest situation and therefore some analytic results can be found [38, 39J. Notice that 
the main qualitative features of the plots in the figures [2] and [3] are independent of the number 
of dimensions. The black curves represent the limiting regimes, which are AdSd+i in the early 
times (bottom curve) and the d + 1 dimensional Schwarzschild black hole at late times (top 
curve). For generic d > 2 the result for AdSd+i is known [7j 



L 



d,reg 



AdS d+1 d — 2 



TV d - "I d ~ 1 
1 \2(d-l)< 



r (2(cZ-l)) 



Unfortunately L^reg for the d + 1 dimensional Schwarzschild black hole is not known. Very few 
analytic results are available for minimal surfaces in four and higher dimensional black holes but 
the curves for Ld ireg have been studied [5H [23] . 

At any intermediate, finite and fixed boundary time t during the black hole formation, we can 
observe from the figures [2] and [3] that L d:ieg (£,t) goes over the Schwarzschild black hole curve 
for small £ and at some point (which depends on t) it leaves from it to adopt a AdSd+i like 
behavior shifted vertically. Indeed, at any finite time the shell is fixed in some region of the bulk 
and we have a AdSd+i geometry inside the shell and a Schwarzschild black hole outside it, with 
the thickness of the transient region is parameterized by a v . For small values of £, the extremal 
surface stays completely outside the shell and far from it feeling therefore only the Schwarzschild 
black hole geometry. As £ increases, the extremal surface begins to enter into the interior part 
of the shell and therefore to feel the AdSd+i geometry inside. This makes L d ^eg deviate from 
the Schwarzschild behavior. When £ is very large, a big part of the extremal surface is inside 
the shell and therefore its length is determined by AdSd+i, explaining the asymptotic behavior 
of the curves in the figures [2] and [3] for large £. The vertical shift in this regime is due to the part 
of the surface which is close to the boundary: being outside the shell, it feels the Schwarzschild 
black hole geometry and it provides a larger contribution to Ld tTe g than the AdS d +i geometry. 
From the plot on the right in the figure [2j we can observe that as the thickness a v decreases, 
Lrf reg reproduces the Schwarzschild black hole result for a larger range of £ and also that the 
vertical shift for large values of I w.r.t. the AdSd+i curve decreases. 

At this point we find it useful to have a look to the shape of the extremal surfaces, which is 
shown in the figure [4] for the case of two disjoint regions in the boundary and will be studied in 
the section [2} The extremal surfaces go backward in the bulk time direction v and around v = 
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they penetrate into the shell probing the AdSd+i geometry inside [33] . 



1.3 Three dimensional case and thin shell limit 



The three dimensional Vaidya geometry (d = 2) is the simplest situation and consequently the 
best one to study in order to get a physical intuition and more analytic results [8j [331 ES ES] 
that could be helpful for the higher dimensional case. Moreover, the CFT on the boundary 
theory is two dimensional and the powerful methods developed for these class of models lead to 
important results [21 El E] which are very helpful to test the holographic techniques. 

For a two dimensional boundary theory, the region A is simply a one dimensional segment 
of length i at fixed boundary time t. The extremal surface we are looking for is given by the 
geodesic connecting the two extrema of this segment and extending in the bulk. The Vaidya 
geometry (1.1) for d = 2 interpolates between AdSs in Poincare coordinates and the BTZ black 



hole of [52]. By employing r 
become respectively [HJ [33] 



and 



1/z as holographic coordinate, the equations (1.14) and (1.16) 



m(v)]{v') 2 + 2rV + r 2 



rv 



<-> i i 
Zr v 



'[0 



.'\2 



1 



. 



(1.21) 



(1.22) 



The first important feature of the three dimensional case is the kind of leading divergence in the 
expansion of the entanglement entropy Sa, which is not power like but logarithmical [21 El [3] - 
Roughly, this could be justified by observing that in the two dimensional boundary theory there 
is no "area law" behavior because the dA is made by two points and it has null measure. In 



three dimensions the regularized area (1.18) coincides with the regularized length of the geodesic 

e = z(£/2 - rj) 



L 



2, reg 




9-V 



dx + 2 log e 



(1.23) 



where e > is the UV cutoff of the holographic direction. The limiting regimes at early and late 
times are respectively AdS^ (bottom black curve in the figure [2]) and BTZ (top black curve) , 
whose regularized lengths read [U [7J 



L 



2, reg 



= 21og£ 

AdS 3 

■ 2ir/VM for I 



L 



2,reg 



2 log 



v smh UJ 



(1.24) 



1. By employing the well known Brown-Henneaux central 



for asymptotically AdSs spaces and the lengths ( 1.24 ) for the geodesies 



where fin = 1/Th 
charge c = 3l/(2G$) 

in ^4^5*3 and in the BTZ black hole, in [61 [7] it was checked that the holographic prescription 
( |1.11[ ) reproduces the expressions for the entanglement entropy of a single interval of length I 
in a two dimensional CFT with central charge c defined on an infinite line at T = and T > 
respectively [21 El E] 



Sa 



T = 



log 



Sa 



T>0 



log 



• V, M 
— smh — 

vre V P 



;i.25) 
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Here (3 = 1/T and the second expression in (1.25) reduces to the first one when /5 
expected. The same happens in (1.24) for M — > 0. 



+oo, as 



Besides AdS% and the BTZ black hole, there is another limiting regime of the Vaidya background 
in three dimensions where the holographic entanglement entropy has been computed analytically 
|38} I39j : the infinitely thin shell limit, defined when the thickness of the shell vanishes a v — > 0. 
In this limit the metric is (1.1) with the mass profile given by a step function m(v) = M 9{v) 
at v = 0, the non vanishing component T vv of the energy- momentum tensor is proportional to 
a delta function 5(v) and the infalling shell represents a shock wave. The resulting geometry is 
given by a planar black brane (a BTZ black hole with a planar horizon) outside the shock wave 
and by AdS% inside the shock wave. The geodesies in both these regimes are known analytically. 
The geodesies entering inside the shell in the full shock wave geometry are piecewise curves 
made by merging the BTZ geodesic outside the shell with the AdS^ geodesic inside. In [381 EH] 
a refraction law has been introduced which tells us how to merge them. Once the merging point 
has been understood, the length of these geodesies is given by the sum of the three pieces: the 
first one going from one extremum of A to the shell in the BTZ geometry, the second one inside 
the shell connecting two merging points of the shell at the same z but different x (AdS^) and a 
third one going from this other point of the shell to the other extremum of A on the boundary, 
again in the BTZ geometry (see the figure [4|. 

Denoting by r H the position of the horizon of the BTZ geometry outside the shell, the regularized 
length of the geodesies entering into the shell can be written as 



L 



2, re; 



2 log 



sinh(r H t) 



where the function s = s(l 
with p defined as follows 



P 



t) 



r H s(£, t) 

sin# G [0, 1] with 9 E [0, 7r/2] can be extracted from 



(1.26) 



2 cos (9 , /2(1 + cos#)p 2 + 2psin6> - cos( 

h log ( — — 

psin# \2(1 + cos9)p 2 — 2psm9 — cos( 



(1.27) 



coth(r H t) + Wcoth (r H t) 



2 cos 6 
1 + cos ( 



(1.28) 



In order to plot -^2,reg as function of £ and of the boundary time t, we need both L2,reg|BTz (1-24) 
and i^reglthin shoii (1-26). Indeed, for small values of i the geodesic is outside the shell and ^2,reg 
is given by .L2,reg|BTz- As I increases, at some point the geodesic enters into the shell and for i 
larger than this value L2, re g is given by L2,reg|thin sheii- From (1.27) one can observe that, for any 



fixed boundary time t, the function I decreases as 9 goes from to 7r/2. Thus, the critical value 
of i from which we have to start using Li re g| t hm sheii is 



=tt/2 



log 



P W=ir/2 



+ 1 



P l0=7r/2 



JL_ / coth(r H t) + 1 
r H ^ \ v coth(r H t) — 1 



2t . 



(1.29) 
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Given a boundary time t, this is the value of I after which the geodesic enters into the shell. In 
the figure [2] the a v = limit we are discussing is represented by the dashed curves. The relation 



(1.29) can be checked on those plots, since the dashed curve characterized by t deviates from 

2i. This is the relation found in 



the BTZ continuous black curve at I = It. This is the relation found in [25] for two dimensional 
CFT models between the duration of the linear increasing of the entanglement entropy after a 
global quench and the size i of the spatial interval A, which lead the authors to suggest the 
quasiparticles picture. 
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Figure 4: Geodesies configuration for the holographic mutual information at the transition point 
for three dimensional Vaidya geometry (d = 2) with M = 1 in (1.9). The total length of the 



geodesies for the connected configuration (red) and the disconnected one (blue) is the same. In 
the upper plots a v = 0.5 while in the bottom ones a v — > (thin shell limit). The boundary time 
is t = 7 (see the intersection of the curves with the horizontal axis in the plots on the right). 
The z axis has been compactified using the arctan function. The green geodesies represent the 
mixed configuration, which is suboptimal. Notice that they do not intersectas it can be clearly 
seen from the plot on the right, top line. 
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2 Holographic mutual information 



In this section we consider the holographic mutual information for three and four dimensional 
Vaidya metrics, employing also the analytic solution for the thin shell limit of the three dimen- 
sional case. We numerically compute its dependence from the size and the distance of the strips 
and also study the transition curves in the configuration space, finding that a time independent 
region exists where the holographic mutual information vanishes at all times. 

When the boundary region A we are interested in is made by two disjoint regions A = A1DA2 
with A\ CiA 2 = 0, the situation becomes more complicated and also more interesting. The most 
important quantity to consider in this case is the mutual information 

I(A U A 2 ) = S Al + S A2 - S AlUA2 (2.1) 

because in this linear combination the leading divergence due to the area law is canceled. The 
quantity I(Ai,A 2 ) measures the classical and quantum correlation between A\ and A 2 . 
Within the context of two dimensional CFT, the Renyi entropies S^ uA = (1— n)^ 1 log Trp AiUA2 
(for interger n 2) have been computed for the free compactified boson [17J and for the Ising 
model |18j . These results have been checked against existing numerical data on spin chains 
[TBI [20] . Unfortunately, the analytic continuation of these quantities for n — > 1 , which is needed 
to find the mutual information through derivation Sa 1 uA 2 = — dnTrp AiUA2 =1 , is still unknown 
for the general expressions, but it has been done for some limiting regimes of the parameters 
like the decompactification regime for the free boson (when the field takes values on the whole 
real line) [17J or when the two intervals are very far apart [18] . The main quantitative lesson 
one learns from these results is that the Renyi entropies, or equivalently the Renyi mutual 
information I^ n \Ai, A 2 ), which is defined in the obvious way by combining the Renyi entropies 



as done in (2.1) for the mutual information, encodes all the data of the CFT. In other words, 
Trp AiUA2 do not contain only the central charge c, like Sa for one interval, but also all the 
conformal dimensions and all the OPE coefficients of the theory. 

From the AdS/CFT point of view, it is natural to study the holographic mutual information, 
which is defined like in (2.1) with Sa given by the holographic formula ( |1.11[ ), as a further test 



of the holographic prescription (1.11) for the entanglement entropy. It is known [14\ I24j that 
the holographic mutual information displays a continuos transition from zero to positive values 
with a discontinuous first derivative which is not observed in the simple CFT models considered 
so far [T71[T8]. This feature is believed to be a large c effect. 

Given the two disconnected regions A\ and A 2 in the boundary, there are three configurations of 
two surfaces extending in the bulk whose boundaries coincide with dA = dA± U dA 2 : the first is 
simply the union of the two surfaces characterizing Sai and Sa 2 (which bounds two disconnected 
volumes and hence it will be referred to as "disconnected" configuration). The second one is the 
"connected" configuration, given by a bridge connecting A\ and A 2 through the bulk (bounding 
a single connected volume in the bulk). The third configuration is composed by the extremal 



13 




Figure 5: Holographic mutual information I(£q,£i,£i) in terms of £q for Vaidya metrics in three 
(plots on the left, infinitely thin shell regime) and four dimensions (plots on the right) in the bulk. 
Different curves are characterized by the boundary time t, whose value increases going from the 
red curves to the blue ones with At = 1 and within the range indicated. The black curves 
correspond to AdS^+i (top curve) and Schwarzschild black hole (bottom curve). The transition 
of the holographic mutual information is continuos with a discontinuous first derivative. 



surface connecting the first extrema of the regions (let us assume for the moment we have strips) 
and the one connecting the second ones. We will denote this case as "mixed" configuration. In 
this mixed configuration, for a static background, the extremal surfaces intersect, but for a 
dynamical background this is not generally true (see the set of green geodesies in the figure |4| . 
Also, for two disjoint circular regions in a 2+1 dimensional boundary, this type configuration 
does not occur at all. A similar mixed configurations of surfaces also occurs in the proof of the 
strong subadditivity for the holographic entanglement entropy [11]. Since they do not always 
intersect for time dependent backgrounds, the proof cannot be extended to the dynamical case 
in a straightforward way. In our case, since Ad^eg is an increasing function of the size of A at 
any fixed time t (see the figures [2] and [3]) we can claim that this mixed configuration is always 



suboptimal with respect to the disconnected one (see also the upper line of the figure 13). Thus, 
the regularized area entering in the holographic computation of Sa 1 uA 2 reads 



Ad,mg - mm (A,reg | connected ' Az.reg | disconn ected ) ' ( 2 ' 2 ) 

Notice that both configurations display the same UV divergence which cancels in the linear 



combination (2.1). In three dimensional backgrounds (d = 2), the lengths of the geodesies are 



involved in (2.2 ). 



It is straightforward to notice that when the disconnected configuration is the minimal one in 



(2.2) the holographic mutual information is zero, while it is positive in the other case. Moreover, 



when A\ and A<i are very far apart from each other the disconnected configuration is minimal. 
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Thus, considering a configuration space which describes all the possible sizes and the relative 
distance between A\ and A%, there must be some region of this space where the "connected" 
configuration is minimal and some other region where, instead, the "disconnected" configuration 
is minimal. The corresponding holographic mutual information is zero and positive respectively. 
The curve in the configuration space which characterizes this transition in the configuration 
space is given by the following equation 

•Ad,reg | connec t e( j "^^i 1 ^? I disconnected ' (^'^) 

This equation is not easy to solve and must be studied case by case. The examples of AdSd+i 
and of the charged black hole in four dimensions have been considered in |14| [24"] . Here we study 
this equation for the Vaidya metrics (1.1) in three and four dimensions. 



t = 1.0* 15 
a, =0.00 
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Figure 6: Holographic mutual information I(£q,£i,£\) as function of l\ at fixed £q for Vaidya 
metrics in three (plots on the left, infinitely thin shell regime) and four dimensions (plots on the 
right) in the bulk. The other parameters are the same ones employed in the figure [5| 



Let us first consider the holographic mutual information I(£q, £1,12) (see (2.1 )) for the Vaidya 



metrics (1.1) in three and four dimensions with the mass profile given by (1.9). This quantity 
depends on many variables and our analysis is mainly numerical. We take equal strips £2 = l\ 
for our plots unless indicated otherwise. 

In the figures [5j [6] we show the dependence of I(£o,£\,£\) from the distance Iq between the 
intervals and the size £\ respectively. The continuos black curves represent the corresponding 
quantities in the limiting regimes of AdSd+i and Schwarzschild black hole (BTZ black hole in 
three dimensions), while the colored ones are characterized by intermediate boundary times 
indicated in the plots. Notice that the qualitative features of the curve do not change with the 
number of dimensions. In general we can clearly observe the transition of the mutual information 
from positive values to zero when £q increases at fixed size £\ (figure [5]) and from zero to positive 
values as £\ increases at fixed separation £q (figure [6]). This transition is continuos with a 
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discontinuous first derivative and unfortunately we do have a clear understanding of it. We 
recall that the holographic mutual information is positive when the connected configuration is 
favored. When we plot a family of curves parameterized by the boundary time t, the common 
feature one observes is that the bigger t is, the larger is the range of variables where the curve 
characterized by t reproduces the corresponding black hole result. Then, for any finite t at some 
point the curve deviates from the black hole behavior and tends asymptotically to the AdS^+i 
behavior, eventually shifted by a constant. 

In three dimensions (d = 2) we can take advantage from the fact that the exact solution is 
known of the thin shell limit case Qy 



(see (1.26)) |38|. I39j and working with the analytic 
solution allows us to extend the range of the variables we can explore. 
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Figure 7: Holographic mutual information I(£o,£i,£i) as function of the boundary time t at 
fixed t\. The different curves are characterized by different values of £q which increase going 
from the top curve to the bottom one. As in the previous figures, we show the d = 2 case on the 
left and the d = 3 case on the right. For a fixed value of £±, varying £q three different behaviors 
are observed. 



As for the dependence on the boundary time t of the holographic mutual information 
I{£o,£i,£{), this is shown in the figure [7] for different values of the separation £q between the 
two strips and a fixed value of their size l\ in three dimensions (plot on the left) and in four 
dimensions (plot on the right). From these plots we observe that, at a fixed of £%, varying 
the separation £q between the strips four different behaviors are observed. For £q very large 
I (£q, £i,£i) is zero at all times. Decreasing £q (i.e. going from the bottom curves to the top ones 
in the figure[7]) we find that I{£q, £\,£\) is zero at t = 0, then it becomes positive for a finite range 
of t and then it vanishes again. Decreasing further £q the holographic mutual information starts 
positive at t = but then it vanishes at some time. Then, for small £q we get that I(£o,£i,£i) 
is positive for any boundary time t, namely the connected configuration is always favored. 
In order to describe these four regimes from another point of view, we find it useful to study 
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Figure 8: Transition point of the holographic mutual information I(£q,£i,£i) for the three 
dimensional Vaidya metric in the configuration space (£i,£q). The black curves represent the 



limiting regimes of AdS% (top curve, given by (2.6) at u = 1) and BTZ (bottom curve, given 
by (2.10) at ui = 1). On the left we plot the transition point for different times (increasing as 
we go from the red curve to the blue one) and thickness a v = 1 (the dashed curve represent the 
thin shell limit a v = 0). On the right, the transition point in the configuration space is plotted 
at a fixed time t = 4 for various values of the thickness a v in (1.9): from a v = (dashed curve) 



to a v = 3 (red curve) . There is a whole region of the configuration space where the holographic 
mutual information is zero for any boundary time. 



the transition curve of I(£q, £i, £2) in the configuration space given by £2, £\ and £q. This means 
to find a family of curves parameterized by t which solves the equation (|2.3|). These transition 



curves are shown in the figures |8j [9] (three dimensional case) and 10 (four dimensional case). 
We set £2 = u>£\ for some finite uj > and then consider the space (£i,£o). In the figure [8] we 
study the three dimensional background for some a v > and uj = 1. In the plot on the right 
the a v dependence is show at fixed t. The dashed curves correspond to the a v — > limit and 
are obtained through the analytic solution of |38|, I39j . The main feature we notice is that it 
is possible to draw a critical curve £o(£\) which is independent of the time t such that for any 
configuration specified by a point above this curve the holographic mutual information is zero 
at all times. This critical curve is above the transition curve of AdS^ and it depends on a v . The 
region below it and above the AdS% transition curve becomes larger as a v becomes smaller, as 
it can be observed from the plot on the right in the figure [8j In the figure [9] we consider this 
transition curves for the analytic solution of the a v — > limit in order to extend the range of 
the configuration space (£±,£0) and to study the dependence on uj, which is equal to 1 as usual 
on the left and equal to 2 on the right. Comparing the two plots we can clearly see the linear 
behavior of the critical curve for large £\ and notice that it depends also on u. 
In the figure 10 we consider the four dimensional background for oj = 1: the qualitative features 
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are the same described above but the configuration space we can explore is much smaller because 
numerically we find it difficult to go further. 




10 20 30 40 10 20 30 40 



t, t, 

Figure 9: Transition point of the holographic mutual information I(£o,£i,u£i) in the configura- 
tion space (£i,£q) for the three dimensional Vaidya metric in the thin shell limit, whose analytic 
solution allows to extend the range of the configuration (see also the figure [8]) . We give the 
curves for different boundary times which increase from the red curve to the blue one. On the 
left we set cj = 1 while on the right oo = 2. The curve above which the holographic mutual 
information vanishes for any boundary time depends on ui. 



2.1 Limiting regimes in three dimensions 



Here we discuss the analytic expressions of the transition curves for the limiting regimes of the 
black hole formation process in three dimensions, namely AdSs for early times and the BTZ 
black hole in the late times (black curves in the figures [5j |6| [8] and [9]) . 

The boundary theory is two dimensional and the two spatial regions A\ and A2 at t = const 
are intervals whose lengths are respectively i\ = X21 and £2 = X43 (we adopt the notation 



1. The separation length is 



X32 and therefore l\ + £2 + £q = X41. 



The transition of the holographic mutual information for AdS^ has been studied in |14j . 
Introducing the harmonic ratio for the four extrema of the two intervals 



Xl2 X34 
X13 X24 



(2.4) 



and employing the first formula in (1.24), the holographic mutual information of two disjoint 
intervals for AdS^ reads Hi 



I(A U A 2 ) 





c 



log 



-1*2 



Uh + h + io) 



log 



x < 1/2 
x > 1/2 



(2.5) 
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Figure 10: Transition point of the holographic mutual information I(£o,ti,ti) for the four 
dimensional Vaidya metric in the configuration space (£i,£q) a t different boundary times which 
increase going from the red curve to the blue one. The black curves represent the limiting regimes 
of AdS^ (top curve) and Schwarzschild black hole (bottom curve). Also in four dimensions there 
is a curve above which the holographic mutual information is zero for any boundary time. 



3Z/(2G (3) 
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The transition point x = 1/2 corresponds to the 

;he logarithm in (2.5) is 
ojI\ with Lj > 0, the equation x = 1/2 becomes a second 



where we recall that c 

solution of hh/i^oih + £2 + £0)} = 1> namely when the argument of the logarithm in (2.5) is 
equal to 1. Parameterizing £2 as £2 

order equation for £q (see (2.4)) with only one positive solution 



to 




4oj 



+ I-IU1 



(2.6) 



(1 + ^) 2 / 

This curve is a line in the plane (£i,£o) passing through the origin whose angular coefficient 
depends on uj. For co = 1 it becomes £q = (\/2 — 1)£\ and this case is employed in the figures [8] 
and [9] (plot on the left). In plot on the right of the figure [9] we set uj = 2. 

The limiting regime al late times is the BTZ black hole. By employing the second formula 



in (1.24), one finds that the equation (2.3) for the transition curve in the configuration space 



given by £2, t\ and £q can be written as follows 

sinh(7rli//?H-) sinh(7rl 2 //3ff) , . 

sinh(7rV/9H)sinh(7r(^i+i 2 + 4)/i8ff) ' 
Introducing u) through £2 = as above and using the addition formulas for the hyperbolic 



functions, the equation (2.7) becomes 

where we have defined 

-— - coth ' 



B(h,uj) + C(£ 1 ,uj)coth 



n£o 



B(£ 1 ,uj) = coth 



+ 1 



C(£t,u)) = coth 



^7 



+ coth 



7TUj£l 



(2i 



(2.9) 



19 



Expressing the hyperbolic functions in their exponential form, the equation (2.8) becomes a 



second order equation in terms of e 2ni °^ H and its positive root provides £q in terms of £\ and 
uj. The result reads 



2tt 



log 



B{tx,u}) + 2 + ,/4[l + B(ei,u)] + C{£i,u>) 2 



(2.10) 



In the figures [8] and [9] this curve is the black one below all the other ones. The curve (2.10) 
passes through the origin (£\,£o) = (0,0) and it always stays below the line (2.6). Moreover, 



the line (2.6) is tangent to (2.10) at the origin and this provides a check of (2.10) because for 
small £\ and finite uj (which implies small £2 as well) the minimal curves remain close to the 
boundary and therefore only the asymptotic geometry of BTZ, which is AdSz, matters. For any 
finite uj > 0, the curve (2.10) tends asymptotically to a horizontal line £q = £q when £\ is large 



for any finite value of uj. In this limit both the disjoint interval are large while the ratio between 
them, being given by uj, is kept fixed. Quantitatively, since for £\ — > +00 we have B(£%,uj) — > 2 
and C{£\,uj) — > 2, the asymptotic value £q reads 



in 



logVl 



7T 



(2.11) 



and it is independent of uj. This means that in the BTZ background, when the separation £q is 
larger than £q, the holographic mutual information is zero for any £\ and £2- 



3 Strong subadditivity and null energy condition 

In this section we explore the relation between the null energy condition for the Vaidya metrics 
and the strong subadditivity condition, which is an important inequality satisfied by the entan- 
glement entropy. We find that a violation of the null energy condition leads to a violation of 
strong subadditivitjj^] 

Consider a quantum system that is partitioned into three or more subsystems, i.e. its Hilbert 
space H can be written as H = ®iHi, and let us denote by Pi lj i 2y „ the reduced density matrix 

obtained by tracing the full density matrix of the system over all Hj with j ^ ii, ii, It can 

be shown, on very general grounds, that the Von Neumann entropy satisfies the subadditivity 
condition 

S( P i) + S(p 2 ) > S( P1 , 2 ) (3.1) 
and also the following two inequalities 

S(pl,2) + S(p 2 ,3) > S( P2 ) + £(^1,2,3) ^ ^ 

S(pi, 2 ) + S( P 2, 3 ) > S( Pl ) + S( P3 ) 
^We remark that this result has been independently obtained also by Robert Callan and Matthew Headrick. 
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Figure 11: Regularized length L reg for three dimensional Vaidya metric when the null energy 
condition is violated: a change of concavity is manifest in both cases. The mass profiles used 
are shown above and they are m(v) = 4f[2 + tanh((u — vo/2)/a v ) — tanh((w + vo/2)/a v )] on the 
right and m(v) = — tanh(-u/a„)] on the left, with M = 1, a v = 0.1, vq = 10. 



which are equivalent and known as strong subadditivity condition (see [HU [TTJ [T2] and the refs 
therein for more detailed discussions). 

If the Hilbert space is partitioned into the product of the Hilbert spaces of local degrees 



of freedom belonging to non intersecting regions of space A\, A%, the inequalities (3.1) and 



(3.2) can be written respectively as 

Sai + Sa 2 ^ Sa 2 uA 2 (3-3) 

and 

Saiua 2 + Sa 2 ua 3 ^ Sa 2 + Sa^a^Az 

(3.4) 

Sa 1 ua 2 + Sa 2 ua 3 Sax + 5*^3 . 

In one dimensional systems (or when the symmetry of the regions considered is such that the 
problem is effectively one-dimensional), for a complete description, it is sufficient to consider 
the entanglement entropy of an interval as a function of its length i, and the two inequalities 
of the strong subadditivity are more conveniently expressed in terms of the function S(£). The 



first inequality in (3.4) states that the function S(£) is concave, and the second that it is non 
decreasing. 



In section [l~T| we mentioned (see ( 1. 10[ ) ) that for the Vaidya metrics ( 1.1 ) the condition m'(v) 



guarantees that the null energy condition is satisfied. By choosing a mass function that does 
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not monotonically increase with v, we can violate the null energy condition and explore the 
consequences of this violation on the entanglement entropy. The results are shown in figure 11 
andE2 
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Figure 12: Location of the transition point of I(£o,£i,£i) when the null energy condition is 
violated: the curve is not monotonically increasing with t\. The mass profiles employed in the 
plots are the same as in figure |TTj 

The curves in figure [TT] are not concave functions of £, but they are still non-decreasing. Therefore 
only the first of the two inequalities is violated. This means that they cannot be equivalent in this 
setting. In order to clarify this apparent contradiction we have to discuss how the equivalence 
between the inequalities is proven, both in quantum mechanics and holographically. 
The two inequalities can be shown to be equivalent by introducing an auxiliary fourth Hilbert 
space H4 such that pi,2,3 = TV4 (t/^I , for a certain pure state [S3]. Then 

S(pi,2,4) = S(p 3 ) and S(p 1A ) = S(p 2 ,3) (3.5) 

and hence, if we write the first inequality for 3 4, 1 o 2 

S(pi, 2 ) + S(p M ) > S( Pl ) + 5(pi, 2j4 ) (3.6) 

and substitute, we get the second, and viceversa. 

If one tries to replicate the argument above in the holographic setting, one encounters a difficulty 
because, although we are guaranteed that it is always possible to find the Hilbert space it 
is not guaranteed that it will be the Hilbert space of the local degrees of freedom of some other 
region of the boundary theory. However, that is the only known kind of partitioning of the 
Hilbert space that allows for a holographic computation. It turns out that, if the bulk manifold 
is homologous to the boundary, the problem is easily solved: can be taken to be the Hilbert 
space of the degrees of freedom of the region = A\ U A 2 U ^3, the complement of A1UA2UA3, 
because it satisfies 

S A S A 2UA3 (3.7) 
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which are the equivalent of (3.6). 



On the other hand, if the manifold contains a black hole, the problem is more complicated, 



because the entropy of the entire system is non-zero, and hence (3.7) do not hold. One has 



to rely on the conjecture that the holographic entanglement entropy is actually describing the 
entanglement entropy of a certain quantum system, and hence on the quantum information 
proof mentioned above. If the dual geometry is sufficiently unphysical, as is the case for when 
the null energy condition is violated, it may not be the holographic description of any quantum 
system, and there is no a priori reason to expect that the two inequalities should be equivalent. 



In terms of the mutual information I(A\, A 2 ) for two disjoint regions (2.1 ) the subaddivity 

I(A 1 ,A 2 )>0 (3.8) 



inequality (3.3) implies that 



while the first strong subaddivity inequality (3.4) can be written as follows 



I(A 1} A 2 UA 3 ) > I(A 1 ,A 2 



(3.9) 



i.e., the mutual information increases as one of the two regions is enlarged while the other one 
is kept fixed. Applying this inequality twice, we one can also conclude that when two equal 
regions are enlarged by the the same quantity, the mutual information increases. 



From the transition curves in figure 12, we can observe that the holographic mutual information 
is not monotonically increasing with £%. This behavior is another manifestation of the violation 
of the first strong subadditivity inequaliy, and has to be contrasted with figure [8] discussed in 
the previous section, where we used the the mass profile (1.9) (figure [I]) and the null energy 
condition is satisfied. 



4 Holographic tripartite information and monogamy 

In this section we consider the holographic tripartite information for three dimensional Vaidya 
metrics and show that the monogamy of the holographic mutual information is violated when 
the null energy condition is not satisfied. 

In addition to the mutual information, another interesting quantity that can be defined from 
the entanglement entropy is the tripartite information 

h(Ai, A 2 , A 3 ) = Sa 1 + Sa 2 + Sa 3 - Sa 1 ua 2 - Sa!ua 3 ~~ Sa 2 ua 3 + Sa 1 vja 2 vja 3 (4.1) 

where A±, A 2 and A 3 are disjoint regions. In contrast with the mutual information, this quan- 
tity is free of divergences even when the regions share their boundary. It is a measure of the 



extensivity of the mutual information; indeed, the definition (4.1) can be written also as 

I 3 (A U A 2 ,A 3 ) = I(A l ,A 2 ) + I{A 1 ,A 3 )-I(A 1 ,A 2 UA 3 ). (4.2) 
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Figure 13: Schematic representations of some mixed configurations occurring in the computation 
of SaiuA 2 (upper line) and 5 , j 4 iU a 2 uA 3 (middle and bottom line). If the regularized area of the 
surface homologous to a single region A is an increasing function of the size of A, then the 
configuration on the left is suboptimal w.r.t. the one on the right (each line on the left should 
be compared with the one on the right having the same color) and thus it does not occur in the 
holographic mutual information or the holographic tripartite information. 

Thus, the mutual information is extensive when 73 = 0, superextensive when 73 < and 
subextensive when 73 > 0. In particular, in either the extensive or the superextensive case, 
namely 

I(A 1} A 2 ) + I(A 1 , A 3 ) < I(A 1 ,A 2 U ,4s) (4.3) 

the mutual information is said to be monogamous. For a generic quantum system, the tripartite 
information can be positive, negative or zero, depending on the choice of the regions. 
Recently it has been shown [12] that for quantum systems with a holographic dual the tripartite 
information is always monogamous. As for the strong subadditivity condition, the proof only 
holds in the case of static dual geometries because in the dynamical backgrounds it is not 
always guaranteed that the surfaces involved in the mixed configurations intersect each other. 
It is therefore interesting to explore the behavior of the holographic tripartite information in the 
simple dynamical backgrounds like the three dimensional Vaidya geometries. 

Among the terms occurring in the definition of the holographic tripartite information, the 
computation of Sa 1 uA 2 uA 3 deserves a short discussion. In one spatial dimension the three regions 
are just intervals, and, according to prescription of [S] for dynamical backgrounds, one has to 
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Figure 14: Schematic geodesies configurations to consider in the computation of SU.1U.A2U.A3- 

find the extremal set of geodesies connecting all the extrema of the intervals. In principle, in 
presence of N intervals one should compare (2N — 1)!! configurations (15 in our case). However, 
since £2,reg(0 at fixed time is an increasing function of £, for iV = 3we are left only with the five 



configurations shown in figure 14, by the argument sketched in the figure 13 Thus, SA 1 uA 2 uA :i 



is given by the minimum among the following quantities 

S(£i) + S{£2) + S{lz) three disconnected volumes 

S(£ 1 +d 2 + £2) + S(d 1 ) + S(£ 3 ) ) 

S(£i) + S(£2 + e?2 + h) + S(d2) } two disconnected volumes 

S(h + d x +£ 2 + d 2 + £ 3 ) + S(d 1 + £ 2 + d 2 ) + S{£ 2 ) J 

S{£\ + d\ + £2 + d 2 + £3) + S(dz) + 5(^2) on e connected volume. 

(4.4) 

Figure 15 displays the time dependence of the tripartite information when the intervals have 
the same size £\ = £ 2 = £3 = 5 and are separated by the same amount d\ = d 2 = £0, with several 
values of £q shown. The behavior is quite complicated and involves different regimes, but the 
quantity is always non positive, even though the geometry is not static. At late times all curves 
go to zero, meaning that, in the thermal state dual to the BTZ black hole, the holographic 
mutual information is extensive. 

Given the results of the section [3] about the relation between the strong subadditivity of the 
holographic entanglement entropy and the null energy condition, we can explore the possible 
relation between the monogamy of the holographic mutual information and the null energy 
condition in the same way, namely by employing mass profiles m(v) which have m'{v) < 
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Figure 15: Time evolution of the tripartite information for the three dimensional (d = 2) Vaidya 
geometry in the thin shell limit. The three intervals have the same size l\ = £2 = £3 = 5 and 
are separated by the same distance d\ = c?2 = £q- The plot shows that the holographic mutual 
information is always monogamous. 



for some range of v. This is relevant because the strong subadditivity and the monogamy are 



indenpendent conditions. In the figure 16 we show the time evolution of the holographic tripartite 

nterval configuratior 
1 at early times (v 



information with the same interval configuration of the figure 15 but with the mass function 
m(v) decreasing from M 



—00) to M = at late times (v — > +00) 
according to the profile shown in the figure 11 (plot on the left). The holographic tripartite 
information becomes positive for certain ranges of t, telling us that a violation of the null energy 
condition leads to a non monogamous holographic mutual information. 



5 Conclusions 

We studied the holographic mutual information for dynamical backgrounds given by the Vaidya 
metrics in three and four dimensions. We found that it is non monotonic as function of the 
boundary time and its behavior depends on the configuration of the two disjoint regions. 
From the transition curves of the holographic mutual information in the configuration space we 
could identify the different behaviors and also find a region in the configuration space where 
the holographic mutual information is zero at all times. Considering the holographic tripartite 
information, we observed that the holographic mutual information is monogamous also for these 
time dependent backgrounds in the ranges of the variables explored. 

By modifying the mass profile occurring in the Vaidya metrics, we showed that the null 
energy condition is a necessary condition both for the strong subadditivity of the holographic 
entanglement entropy and for the monogamy of the holographic mutual information. A deeper 
understanding of the relation between the null energy condition and the inequalities satisfied by 
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Figure 16: Time evolution of the tripartite information for the three dimensional (d = 2) 
Vaidya geometry when the null energy condition is violated: the mass function decreases from 
m(— oo) = 1 to m(+oo) = 0, according to the profile shown in the figure 11. Since I3 becomes 



positive for certain range of t, the monogamy of the holographic mutual information is violated. 



the quantities defined from the holographic entanglement entropy is needed. 
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